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SUMMARY 


This paper describes the initial results of an effort to develop partial 
coherence technigues for interior noise source/path determination in the highly 
coherent environment of propeller-driven general aviation aircraft. Examples 
illustrate the effects of measurement interference and the use of a two-channel, 
real-time analyzer for the analysis. The paper includes a summary of the 
computational techniques and illustrates their application to a two input, 
single output system with coherence between the inputs. Errors introduced into 
the calculations by the method used for data analysis are discussed. The 
results illustrate the importance of using a simultaneous time base for the 
data reduction and indicate the type of errors that can be encountered by 
failure to observe this requirement. 





INTRODUCTION 


Effective procedures for controlling interior noise require identification of 
the noise sources and the noise transmission paths. Current methods for 
accomplishing this are cumbersome and time consuming. The difficulties 
encountered are particularly evident in the highly coherent environment of a 
propeller-driven general aviation aircraft, A new technique for efficient, 
reliable determination of noise sources and paths along with the capability to 
rank order their importance is needed. 

Recent developments in computational procedures have led to increased interest 
in partial coherence analyses for source/path determination. The theoretical 
aspects of this method as developed by Dodds and Robson (ref. 1) and further 
improved by Bendat (refs, 2, 3, 4) are particularly appropriate for tlv's 
investigation. Applications of these methods to diesel engines (refs. 6), 
a punch press (ref. 7), and a light aircraft (ref. 8) have been reported in the 
literature. However, these previous authors reported only partial success in 
identifying sources. Improvements of the approach are felt to be needed before 
the method can be successfully applied to interior noise problems encountered 
in propeller-driven general aviation aircraft. 

The purpose of this paper is to describe the latest results of an ongoing effort 
to develop partial coherence techniques for interior noise source/path determi- 
nation in the highly coherent environment of propeller-driven general aviation 
aircraft. The paper includes a summary of the theoretical method as developed 
by Bendat (ref. 2) and illustrates the application to a two input, single output 
system with coherence between the inputs. The augmentation of the calculations 
on a digital computer interfaced with a two-channel real-time analyzer is 
discussed. The results presented indicate possible sources of error in the 
computations and suggest procedures for avoiding these errors. 


DESCRIPTION OF ANALYSIS 


A schematic indicating the analytical model of a physical system used in partial 
coherence analyses is shown in figure 1, The physical inputs consists of a 
number of time histories which have been directly measured on the physical 
system under consideration and are assumed to represent the various sources of 
interior noise. There may be various degrees of coherence between these input 
records as well as between each of them and the output record which has been 
directly measured also and is assumed to represent the interior noise environ- 
ment (receiver). It is assumed that none of these coherences are identically 
one. If this occurs, the records contain redundant (or unnecessary) information 
and some of the records should be eliminated from the analysis (ref. 2). After 
the redundant records are eliminated, the remaining input records are ordered. 
Although the selection of the order of these records is largely arbitrary, one 
procedure is to choose the record with the highest coherence between the input 
and output as the first ordered input; the record with the next highest 
coherence is chosen as the second ordered input, and so forth. 
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Once the ordered inputs have been obtained, the rest of the analysis is usually 
carried out in the frequency domain. Conceptual ly> however, the process is 
equivalent to obtaining conditioned inputs in the time domain. The ti'^st 
conditioned input is identical to the first ordered input, The second 
conditioned input, X 2 .-|(t), is the second ordered input with the effects of 

X](t) removed. The third conditioned input, X3,-|2(t)> is the third ordered 

input with the effects of xi(t) and X£(t) removed, and 
inputs are then mutually uncoherent. Each conditioned 
to be the input signal to an ordered transfer function 
particular conditioned input record to the output. The ordered transfer 
functions are not unique, but depend upon the particular order used the 
ordered inputs. Equations relating the physical transfer functions to the 
ordered transfer functions may be found in reference 3. The relative effect 
each ordered input on the output is given by an equation of the form 
Syy = 1h12Sxx (ref. 4). 


on. The conditioned 

input record is assumed 
which relates that 


of 


As previously stated, the analysis is usually carried out in the frequency 
domain. The first step is the calculation of the auto-spectra and cross- 
spectra for all possible combinations of the ordered inputs and the output. 
After these spectral functions have been calculated, the effects of the inpu 
are removed from the spectra. The formulas for removing the effect of 
input k from the auto-spectra and cross-spectra are shown in figure 2. As 
indicated in the figure, the equation for (auto-spectrum of ordered input 


with the effect of input k removed) has two forms 
which explicitly involves the coherence function 


The form of this equation 


y 2 , clearly demonstrates 
ik 


that > 0 since < i • However, the other form of this equation 

involves multiplications, divisions, and subtractions of three different 
spectra and it is entirely possible that small, statistically insignificant 
differences in the estimates of these spectra could result in the computation 
of negative numbers. The formula for (cross-spectrum between ordered 

inputs i and j with the effect of input k removed) involves four different 
spectra and is also subject to the kinds of errors just mentioned. The 
formulas for removing the effects of more than one input are similar and may be 
found in the literature, e.g. , reference 2. 


Also shown in figure 2 is the partial coherence function yfy.^ (pai^tial 

coherence between ordered input i and the output with the effect of input k 
removed). Note that partial coherence is simply ordinary coherence computed 
using conditioned spectra and may be used to rank order the importance of the 
inputs. That is, it follows from the formula for the output due to ordered 
input i with the effect of input k removed (fig. 2) that Y^y,|, < ^jy-k 

implies ordered input j produces a greater part of the output signal than 
ordered input i (with input k removed). 
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RESULTS AND DISCUSSION 


Source Measurement Interference 

The problem of source measurement interference provides a typical application 
for partial coherence techniques. The analytical model for this problem is 
shown in figure 3, The two sources of noise are assumed to be uncoherent. 

Two microphones, M] and M 2 » a'f'e used to measure the noise sources but because 
of measurement interference, the quantities actually measured are 
Mi: Sii + k]S22 and M 2 : S 22 + k 2 Sn, where the parameters k] , k 2 determine the 
coherence between the input measurements. For the numerical example presented 
herein, the coherence between the input measurements is y2^ = 0.7, An 

analytical transfer function is specified between each of the source inputs 
and the output. The auto-spectra shown in figure 3, along with the various 
cross-spectra, were used in the computational algorithms to compute the partial 
coherence functions between M] , M 2 , and the output, and estimate the previously 

specified transfer functions. Figure 4 shows plots of two of the coherence 
functions for this test case. The ordinary coherence function between input 1 
and the output is nearly equal to 1 for frequencies up to 500 Hz and drops to 
a value of about 0.8 for frequencies between 700 Hz and 1000 Hz. The partial 
coherence function between input 1 and the output with the effect of input 2 
removed is identically 1. This value is correct since, after the effects 
of input 2 are removed, the model consists of a single input, single output, 
linear system. 

Figure 5 shows a comparison of the exact and estimated transfer functions for 
the test case. The first resonant peak of Hiy at 400 Hz is determined 

reasonably well by the estimated transfer function. The frequency of the 
second resonant peak of H-j^ at 700 Hz is also determined accurately but the 

amplitude of the transfer function is somewhat overestimated at this frequency. 
For H 2 y, the resonant peak at 400 Hz is also determined quite well by the 

estimated transfer function. However, at 700 Hz the estimated function 
indicates a low amplitude resonant peak which is not present in the exact 
expression for H 2 y. The presence of this peak in the estimated transfer 

function is due partly to the coherence which exists between the inputs, M] 

and M 2 , and partly to the difference between ordered transfer functions and 
physical transfer functions. This result indicates that care must be exercised 
in inferring the existence of resonant responses from a partial coherence 
analysis using experimental data for which measurement interference may exist, 
and for which the actual resonances are unknown. Although the amplitudes of 
the estimated transfer functions are somewhat in error, these transfer functions 
do correctly indicate the relative importance of the two inputs in this example. 
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Experiment Using Computer^'Analyzer Interface 

In this test case, the use of a hard-wired, two-channel analyzer was investigated, 
A set of experimental data was obtained by simultaneously ^ 

of a white noise generator on three channels of a tape recorder. These data 
Sere Snsidered as a two input, single output system with transfer f^^nctions of 
unity between each of the input signals and the output. The data were analyzed 
on a^desk top computer which was interfaced with a two-channel analyzer using 
the IEEE bus. This interface was straightforward and presented no particular 
oroblems However, the use of a two-channel analyzer does have implications on 
tJlfsSEs^quenrS analysis. This system dictated that only three spectra 
could be obtained simultaneously: two auto-spectra and che * . 

Thus, several passes of the data were required to obtain the necessary spectra. 

Since it was impossible with the equi^ent a^tuafLearf 

at precisely the same spot on the analog tape f^ery time, the ^‘^tual sp 
used in the analysis were obtained from slightly different time segments. This 
w^found to adversely affect the accuracy of the results specifically in that 
values of coherence* much greater than 1 were computed. As a measure of this 
accuracy loss, the ordinary coherence functions have been recomputed using the 
speara stored in the compiter and compared with the coherence functions from 
the two-channel analyzer which were obtained as the data were processed. It 
is felt that any additional results obtained by futher processing of this data 
will be no more accurate than these recomputed coherence functions. 

Figure 6 shows several examples of these recomputed ordinary coherence 
functions. As indicated in figure 6(a), which was recomputed using three 

spectra obtained simultaneously, displays the characteristics expected for 
these coherence functions. Its value is very nearly 1 over most of the 
frequency range. All of the coherence functions for this case should appear 
vir?uallj Identical to the coherence function rfj ^'^ure, if 

the computations are accurate. 

The coherence function y^ shown in figure 6(a), which was recomputed with 

only two spectra obtained simultaneously, indicates errors of about ±10 P^'f'cent 
(compared with y^^). Some of the values are greater than 1 and could lead to 

the computation of negative values for the 

figure 2.) The coherence function y^^ recomputed with none of the spectra 

from simultaneous time segments, shows large deviations (±30 percent) from the 
true value. Errors of this magnitude obviate the usefulness of coherence 
analyses. 

Figures 6(b) and 6(c) show the effects of increasing the number of averages 
(figure 6(b)) and increasing the bandwidth of the analysis (figure 6(c)) for 
a recomputed coherence function with two spectra obtained simultaneously. 

As these figures indicate, increasing either the number of averages or the 
bandwidth of the analysis does have a smoothing effect on the recomputed 
coherence functions. However, the amount of data available for analysis or 
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the frequency resolution required to separate different sources may easily make 
either of these approaches unfeasible. 

An error analysis was conducted to determine the source of the loss of accuracy. 
Standard error formulas were used to obtain 95 percent confidence intervals 
for the individual spectra (ref, 4), The worst combinations of the numbers so 
obtained were then used in the equation for calculating coherence functions to 
obtain the actual errors. The results of this analysis are shown in figure 7, 

The solid curve in figure 7 shows the percentage of the error (which may be 
plus or minus) versus the number of averages for data with two spectra obtained 
simultaneously. The dashed curve shows the error for data with none of the 
spectra obtained simultaneously. Also shown on figure 7 are experimental 
results for two spectra obtained simultaneously (circles) and none of the spectra 
obtained simultaneously (squares). As the figure shows, the experimental results 
are in good agreement with this error analysis. 

Although the data analysis can be accomplished using a two-channel analyzer with 
no fewer than two spectra obtained simultaneously, such an approach would 
require reanalyzing the same data several times. This fact, together with the 
large number of averages required to maintain reasonably small errors (say 
10 percent), suggests that alternative methods of analysis, such as simultaneous 
digitizing of all data, would be preferred. In addition, the results of this 
study indicate that the calculations which must be performed in a partial 
coherence analysis are quite sensitive to such parameters as the number of 
averages and the bandwidth of the analysis. 


CONCLUDING REMARKS 


This paper has described the initial results of an effort to develop partial 
coherence techniques for interior noise source/path determination in the highly 
coherent environment of propeller-driven general aviation aircraft. Examples 
are shown to illustrate the effects of measurement interference and the use of 
a two-channel, real-time analyzer to obtain the initial spectra required for the 
analysis. Use of the two-channel analyzer system illustrates the importance of 
obtaining all of the necessary spectra from simultaneous time segments. Failure 
to observe this requirement introduces errors of up to ±30 percent into the 
computations. Such errors may result in computed coherences much larger than 1. 
Results were presented showing that these errors could be reduced by using a 
large number of averages, up to 2000, in obtaining the spectral estimates. 
Because such a large number of averages is difficult or impossible to obtain in 
practical experimental situations, alternative methods of analysis, such as the 
simultaneous digitizing of all data, are preferred. 

An analytical example of the measurement interference problem has also been 
presented. The results indicate that the estimated transfer functions can be 
used to determine the relative importance of the noise source inputs. 
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Figure 1.- Schematic of partial coherence analytical model for multiple input system. 
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Figure 2.- Calculation of conditioned spectra and partial coherence functions. 









Figure 4.- Coherence functions for source 
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Figure 5.- Transfer functions for source measurement interference problem. 
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Figure 6.- Examples of accuracy loss in recomputed coherence due to nonsimultaneous spectra. 
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Figure 6.- Concluded. 
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Figure 7.- Error due to nonsimul taneous spectra. 
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